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Abstract 



We classify the niatrix product states having only spin-flip and parity syni- 
5 I metries, which can be constructed from two dimensional auxiliary matrices. We 

show that there are three distinct classes of such states and in each case, we 
^^1 determine the parent Hamiltonian and the points of possible quantum phase 

transitions. For two of the models, the interactions are three-body and for one 

the interaction is two-body. 
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1 Introduction 

The problem of determining the ground state of a given many-body Hamiltonian, 
is an important problem in condensed matter and mathematical physics. There is 
already a rich literature on this subject, which dates back to the work of Hans Bethe 
on the Heisenberg spin chain and has continued since then with the works of many 
other people including Yang, Baxter, and Lieb, to name only a few. In particular 
for spin systems, the exponential increase in the dimension of Hilbert space of such 
a system, as the number of particles rise, turns this problem into a computationally 
formidable one, beyond the capability of any classical computer. In fact it has now 
been established that finding the ground state of a given many-body Hamiltonian is 
the analog of NP-complete problems for quantum computers [1] . The lesson that we 
learn from all this is that it is highly improbable that we be able to find a generic 
system with exactly known ground state. Nevertheless, there are systems with ex- 
actly known ground states and even if such systems are not exactly what we have in 
nature or in the laboratories, they may be good approximations to real systems, or 
at least may teach us useful and important concepts and methods for studying more 
realistic systems. 

One of the methods, developed in recent years, for investigating this problem is 
the Matrix Product State or Finitely Correlated State [21 O Hj formalism. It is also 
called Optimal Ground State formalism in some references O |6l [7] . The main theme 
is that one starts with a state with prescribed symmetries and properties, and then 
construct the family of Hamiltonians for which this state is an exact ground state. 
It is obvious that for any given state \ip), the equation H\'iIj) = has always many 
solutions for the unknown H, since the number of equations is much less than the 
number of unknowns. However the problem becomes interesting and quite non-trivial 
when we put physical constraints on the Hamiltonian. That is we demand that i) 
H be positive, so that |^) is actually the ground state and not an ordinary eigen- 
state, ii) that it be a sum of local terms, i.e. H = J2k hk,k+u where hk,k+i acts on a 
block of / -|- 1 spins, and iii) that both the state and the Hamiltonian have some rea- 
sonable physical symmetries, like parity, spin-flip, and at times rotational symmetries. 

In the past few years, a lot of interest has been attracted to the subject of matrix 
product states [H |9l [lOl [TTl [121 US] , specially after the emergence of the field of quan- 
tum information [141 [151 IISl [1^ 118] . The reason is the complementary role that the 
fields of condensed matter physics and quantum information play in investigation of 
many body systems. On the one hand quantum information starts with properties 
of states, while condensed matter physics, starts from the properties of the Hamilto- 
nian which embodies the interactions and energy of the system. The matrix product 
formalism is one of the subjects which lies at the borderline of these two subjects. 

As is well known, in this formalism, one starts from proposed states whose ex- 
pansion coefficients are the trace of product of given matrices. While for numerical 
investigations, i.e. the density matrix renormalization group (DMRG), one usually 



starts from large dimensional matrices, to simulate ground states of given Hamilto- 
nians, in the approach which is used for finding exactly solvable models, one starts 
from low dimensional matrices and finds Hamiltonians for which these states are exact 
ground states. This is the approach which has been used in our works and in many 
other works in the past few years [SlElIinillllliaiElIiaEniEIlESEaiMlES]- In 
this article we follow this approach and classify all the matrix product states which 
can be constructed from two dimensional matrices. We restrict ourselves to states 
which allow one or another of the spin-flip or parity symmetries and find that there 
are three classes of such matrix product states. We will study these states and find 
the parent Hamiltonians and also the points or lines in the space of control parame- 
ters where a MPS-quantum phase transition ^26j (MPS-QPT) may occur. This is a 
term, introduced in [27] to differentiate these kinds of QPT's (characterized by any 
discontinuity in any macroscopic quantity) from the conventional QPT,s in which a 
non-analyticity in the ground state energy typically occurs. 

The structure of this paper is as follows: To make the article self-contained, in the 
next section we briefly introduce the basic elements of the formalism. In section ([3]) 
we discuss the symmetry properties of MPS and classify the spin 1/2 states with two 
dimensional matrices. We show that there are three classes, denoted by a model A, 
model B and model C and studied in subsequent sections. We end the paper with a 
discussion. 

2 A brief introduction to matrix product states 

First let us review the basics of matrix product states. Consider a homogeneous ring 
of N sites, where each site describes a d— level state. The Hilbert space of each site 
is spanned by the basis vectors |i), i = 0, • • • d — 1. A state 

I*) = Yl V'Jii2-iivKl'^2,---,^Ar) (1) 

is called a matrix product state if there exists D dimensional complex matrices Ai E 
C7^x^, i = 0---d-l such that 



where Z is a normalization constant given by 

Z = tr{E^) (3) 

and 

E:=Y,A*®U (4) 

i=0 

Here we are restricting ourselves to translationally invariant states, by taking the 
matrices to be site- independent. By defining the vector valued matrix 

d 

A = Y.Ai\i), (5) 



one can write the MPS in a more concise way as 

\iIj) = tr{A(g>A(g>---A), (6) 

where we use the convention 

tr{A ^ A) := tr{AiAj)\i) ^ \j). 

It is important to note that the MPS representation ([2]) is not unique and a 
transformation such as 

Ai -^ fiUAiU-^ (7) 

where U is an invertible matrix, and /x is a constant, leaves the state invariant. The 
simple structure of the MPS allows also an easy calculation of correlation functions. 
Let O be any local operator acting on a single site. Then we can obtain the one-point 
function on site k of the chain (^|0(A;)|'I') as follows: 

„io(.)i*, . "'";;y , (8) 

where 

d-i 

Eo:= J2ii\0\J)A*'^A,. (9) 

i,j=0 

In the thermodynamic limit A'^ — > oo, equation ([H]) gives 

(^|0|*) = ^^'""^l^ol^""'^) ^ (10) 

■^max 

where we have used the translation invariance of the model and Xmax is the eigenvalue 
of E with the largest absolute value and \Xmax) and {Xmax\ are the right and left 
eigenvectors corresponding to this eigenvalue, normalized such that {Xmax\Xmax) = 1- 
Here we are assuming that the largest eigenvalue of E is non-degenerate. In case X^ax 
is degenerate with degree equal to g, then Eq. (fTOj) will be modified to 

('^lOl'^) — ■^»=l \Xmax,i\Eo \ Xmax,i) /-. -, \ 

'^max 

The n-point functions can be obtained in a similar way. For example, the two-point 
function (^|0(A;)0(Q|^) can be obtained as 

moimm - '^'^t^ggj"^"' (12) 

where Eo{k) := E EqE^ . Note that this is a formal notation which allows us 
to write the n-point functions in a uniform way, it does not require that E be an 
invertible matrix. In the thermodynamic limit the two point function turns out to be 

(^|0(l)0(r)|M/) = J—Y^Xl-^{X^,,\Eo\Xi){X,\Eo\X^ax)- (13) 



For large distances r ^ 1, this formula reduces to 

xr-2 

(^|0(l)0(r)|^) - {^\0\^f = -^{X^,,\Eo\Xi){Xi\Eo\Xma.), (14) 

where Ai is the second largest eigenvalue of E for which the matrix element ( Ai | Eq \ Xmaa 
is non-vanishing and we have assumed that the eigenvectors of E have been normal- 
ized, i.e. (Aj|Aj) = 5ij. Thus the correlation length is given by 

( = ^^- (15) 

Any level crossing in the largest eigenvalue of the matrix E signals a possible MPS- 
QPT . Here we are using the term quantum phase transition in a broader sense than 
usual [27], that is, we call any discontinuity in any macroscopic quantity a quantum 
phase transition, even if the ground state energy itself is a continuous function of the 
coupling constants. Also, due to (llSp . any level crossing in the second largest eigen- 
value of E implies the correlation length of the system has undergone a discontinuous 
change. 

2.1 The Hamiltonian 

Given a matrix product state, the reduced density matrix of k consecutive sites is 
given by 

triiA* ■■■A* ® A^, ■ ■ ■ A,- )E^-^) 
Pn-ikjvjk - tr(E^) ' 

The null space of this reduced density matrix includes the solutions of the following 
system of equations 

d-l 

E cj,..,,^,v^i. =0. (17) 

Given that the matrices Ai are of size DxD, there are D^ equations with d unknowns. 
Since there can be at most D^ independent equations, there are at least d^ — D^ 
solutions for this system of equations. Thus for the density matrix of k sites to have 
a null space it is sufficient that the following inequality holds 

d'' > D'^. (18) 

Let the null space of the reduced density matrix be spanned by the orthogonal vectors 
Icq,), {a = 1,- ■ ■ s,> d^ — D^). Then we can construct the local hamiltonian acting 
on k consecutive sites as 

s 

h := ^fJ.a\ea){ea\, (19) 

j where /ia's are positive constants. These parameters together with the parameters 
of the vectors \ei) inherited from those of the original matrices Ai, determine the total 
number of coupling constants of the Hamiltonian. If we call the embedding of this 



local Hamiltonian into the sites I to I + k hy hi^i^k then the full Hamiltonian on the 
chain is written as 

N 

H = Y,hi,i+k- (20) 

1=1 

The state |^) is then a ground state of this hamiltonian with vanishing energy. The 
reason is as follows: 

N 
{^\H\^) = tr{H\^)m = J2 tr{hi,l+kPl,l+k) = 0, (21) 

1=1 

where pi^k+i is the reduced density matrix of sites I to l + k and in the last line we have 
used the fact that h is constructed from the null eigenvectors of p for k consecutive 
sites. Given that H is a positive operator, this proves the assertion. 

In view of the above introduction, we have a clear recipe for constructing matrix 
product states and a family of parent Hamiltonians. First one chooses the matrices 
throwing away all spurious degrees of freedom by transformations ([7]) and reducing 
further the degrees of freedom by imposing symmetries. In this way one ends with 
a reasonable set of matrix product states, which hopefully may have applications in 
description of real physical systems. Imposing a continuous symmetry, like rotation 
around an axis, restricts the matrices considerably [8l |9l [10]. In this article we re- 
strict ourselves to discrete symmetries only which allow a larger variety of models 
to be constructed. For two dimensional auxiliary matrices, this is a simple tractable 
problem, which we do in this article. For larger matrices, the problem is not so simple 
and we defer it to another work. 



3 The classification of 2 dimensional matrices for matrix 
product states of spin 1/2 chains 

We now classify all the two dimensional matrices which can be used for constructing 
spin 1/2 matrix product states. We restrict ourselves to the case where these states 
have spin- flip and left-right symmetries. 

3.1 Symmetries 

Consider now a local symmetry operator R acting on a site as R\i) = Rji\j) where 
summation convention is being used. R is a d dimensional unitary representation of 
the symmetry. A global symmetry operator TZ := R'^ will then change this state to 
another matrix product state 

^i,,,...i^ -^ V := ir« A^^ • • • A', J, (22) 

where 



A'i ■= RijAj. (23) 



A sufficient but not necessary condition for the state |^) to be invariant under this 
symmetry is that there exist an operator U{R) such that 

R^jAj = U{R)AiU~^{R). (24) 

Thus R and U{R) are two unitary representations of the symmetry, respectively of 
dimensions d and D. Equation (|24p wih be our guiding lines in defining states with 
prescribed symmetries. Spin-flip symmetry means that 

where i = 1 — i. For a matrix product state, this requires that there be a matrix like 
X, such that 

XAoX~^ = eAi, XAiX-^ = cAq, (26) 

where e = ±1. Similarly left-right symmetry means that 

For a matrix product state, this means that there be a matrix Q such that 

QAon-^ = cjA^, ftAiQ-^ = aAj, (28) 

where the superscript T stands for the transpose and o" = ±1. These conditions 
are general irrespective of the dimension of matrices. For two dimensional matrices 
however, if we take the trace and determinants of both sides of equations (I26p and 
(j28p . and comparing them, we find that 



tr{Ao) = etr{Ai), (29) 

and 

det(^o) = det{Ai). (30) 

The important point is that for two dimensional matrices, the trace and determi- 
nant are the only invariants under similarity transformations, and hence these two 
equations allow us to classify all the matrices ^o ^-nd Ai which can be used for con- 
struction of spin 1/2 matrix product states. We will use the freedom ([7]) and also the 
above two conditions to show that there are three distinct classes of matrix pairs and 
corresponding matrix product states. In the next three sections, we introduce the 
matrix pairs and study the properties of the matrices obtained from them. For ease 
of distinction, we use different notations for the matrix pairs in each section, namely 
we denote the matrix pairs by Ai, Bi and Cj in the following sections. 

4 Model A 

If one of the matrices say Aq is diagonalizable, we can use freedom in re-scaling 
Aq — > fiAQ to put it in the form 



^o = r:^ , 1, (31) 




where g is a, free parameter. Now take Ai as an arbitrary matrix of the form Ai 
I I . If 6c 7^ 0, then we can use the transformation ([7]) with f7 = /r 



and a further re-definition of Vbc — > b to put it in the form 

M = (l ']■ (32) 



For this type of matrices, the MPS is automatically left-right symmetric, with fi = /. 
From equations (j29p and ()30p . we find the following constraints on the parameters, 

I - g^ = (^ad - b^) 2e = {a + d). (33) 

To solve the second equation, we put 

a = e + u, b = e — u, (34) 

which turns the first equation into 

which can be solved by the parametrization u = g cos 6 and b = g sin 6. Therefore the 
final form of the matrices will be as follows: 

Ao=('1,' ' ], A, = (^ + ^.7^ ^^"^^„V (36) 

\ 1 — g I \ gsmO €. — g cos 6 j 

The matrices satisfy the symmetry constraints (j26p and ()28p with Vi = I and 

\ 1 — e cos 1 + e cos 6 j 

We now restrict ourselves to the case e = 1 (We do this also for other models B 
and C). To study the properties of the corresponding MPS, we should determine the 
eigenvalues of the transition matrix Ea '■= Aq ® Aq + Ai ® Ai. For general values of 
the parameters, the analytical form of these eigenvalues are complicated. However we 
can gain insight by looking at them for for generic values of the parameter g. Figure 
^ shows the eigenvalues as a function of 6 for two values of the parameter g. The 
same pattern repeats for other values of g. We first see that there is no level-crossing 
in the largest eigenvalue and hence no MPS quantum phase transition in this model. 
However for every value of g, two points are important. The point 9 = 11 (or — H), is 
the point where the largest and the next to largest eigenvalues become equal. This 
is a point where according to (llSp . the correlation length becomes infinite. It is seen 
from ()36p that at this point, the matrices become diagonal Aq = diagonal{l + g, 1 — g) 
and Ai = diagonal{l — g, l+g) and according to ([6]), the un-normalized state becomes 
the sum of two product states, namely 

IV^) = |</)+)«^ + |(/)_)®^, (38) 








Figure 1: (Color Online) The eigenvalues of the transition matrix Ea for e = 1 and a 
generic value of g (i.e. g=l/2) as a function of 9. 

where 

!</>+) = (1 + g)\0) + (1 - g)\l), !</._) = (1 - gm + (1 + g)\l). (39) 

Let us now find the parent Hamiltonian for a fixed value of the parameter g, say 
g = 1. To this end, we have to see for which value of k (the range of interaction), the 
system of equations (|17p have a non-trivial solution. It is seen that the smallest k for 
which there is such a solution is A; = 3. 

The solution space of this system of equations turn out to be spanned by the 
following vectors: 



\ef) 

\et) 
let) 



1 + cos(6l) 



|000) + |101), 



|001) 
1 100) 



|011) 
|110) 



1 + cos(9) , 
|010) ^|111) 



(40) 

(41) 
(42) 

(43) 



The symmetries of the parent Hamiltonian now show itself in the form that the 
above states, which transform into each other under the action of these symmetries. 
To have a Hamiltonian which respects these symmetries, we construct it as follows: 



hA = J{\et){et\ + |ef )(ef I) + K{\ei){e^\ + |e^)(e^|). (44) 

The final form of the the full Hamiltonian in terms of Pauli matrices, after neglecting 
additive and multiplicative constants becomes 



N 

Ha = Y. ^i^^f+i + 'J^'^K+2 - uJ(Tt(Tt+2 + uJcy\ay+2 

(45) 

(46) 
(47) 

This is a three-body Hamiltonian with two free coupling constants. We will say more 
about this in the discussion. 

5 Model B 

In accordance with our proposed notation, we denote the matrices in this case by Bq 
and Bi. In this case the matrices are the same as in the previous case, except that 
one of the parameters, say h is zero. There is no transformation which can put Bi 
into symmetric form, and we have 



K r ^ z r z 




1+^2°'' and 








---'^'V- 


K 




^1 = : . • m 



In this case no similarity transformation can put the matrix B\ into symmetric form, 
without destroying the diagonal form of Bq and hence the MPS will not be parity 
invariant or left-right symmetric. From Eqs. (j29p and ()30p . we find that 



\-g^ = ad 2e = a + d, (49) 

which lead to the following final parametrization for the matrices, where for definite- 
ness we will show the matrices by a different letter 



Bo=ri' .\i ^1= : . . • (50) 





In this case again we have two free parameters in the MPS state, namely c and g. 
Moreover the matrices satisfy the spin flip symmetry condition ()26p with 




X = ^ ^ +^ ^ I (51) 



while they do not have any symmetry under parity (i.e. there is no matrix 0, satisfy- 
ing (EiD). 

The eigenvalues of the transition matrix Eb = Bq <^ Bq + Bi Bi, (for the case 
e = 1) now become 

Af = 2(l+5)^ Af = 2(l-g)2, Xf, = 2{l-g^), (52) 

independent of the value of c. Figure ([5]), shows these eigenvalues as a function 
of the parameter g. From this figure, a few features can be recognized. First we 
note that at g = 0, there is a crossover between the largest and the second-largest 
eigenvalues. This points to a possible MPS-quantum phase transition at this point. 
Furthermore at g = ±1, there is a discontinuity in the derivative of the second largest 
eigenvalue which points to a discontinuity in the derivative of the correlation length. 
In view of the general and rather broad definition of MPS quantum phase transition, 
as the appearance of any discontinuity of a macroscopic observable, these points are 
also points of MPS-QPT's. Moreover since the eigenvalues do not depend on the 
parameter c, it appears that the above points are really lines in the space of coupling 
constants c and g. 

From ()36|) . it is clear that at g = 0, the un-normalized MPS turns into the following 
state 

l^i?) = lx+r^ + IX-)^^ (53) 

where 

|X±> = (1±5)|0) + |1). (54) 

Finally we come to the parent Hamiltonian. For this model we find that the 
smallest value of k for which the system of equations (J17l) has a non-trivial solution 
is k = 2 and hence we can have a two-local parent Hamiltonian. The solution space 
of the system of equations (J17p is spanned by the vectors 

|ef) = -^(1 + 5)|00) + |01) + 1(5-1)111) (55) 

|ef) = -^(1 + 5)|11) + |10) + ^(5-1)100). (56) 

Interestingly we note that the above vectors transform into each other under spin- 
flip, but they do not have any transformation property under parity which is to 
be expected, since the original matrices had only spin-flip symmetry. The parent 
Hamiltonian which is symmetric under spin flip will be 
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Figure 2: (Color Online) The absolute values of the eigenvalues of the transition 
matrix Eb for e = 1, as a function of g. 



B\/R\ 



hB = J{\ef){ef\ + \e^){e^\) 



(57) 



and the full Hamiltonian will be (after neglecting additive and multiplicative constants 
and collecting all the various terms) 



TV 



H^ = J2il-ff')(afat^,-aya, 



i+l. 






(58) 



This is the Heisenberg XYZ system with specific couplings, that is we have found 
exact solution on a submanifold in the space of couplings. Note that since the Hamil- 
tonian does not depend on the paramter c, while the MPS does, this means that 
there is a large degeneracy in the ground state. Expansion of the MPS in terms of 
the parameter c, i.e. |^(c, 5)) = J2n c"|V'n(5)), will yield the multitude of degenerate 
ground states \ipn{g)) as in [28] . 

6 Model C 



Denoting the matrices by Co and Ci, this is the only remaining case, where Cq is not 
diagonalizable and can be put only in the Jordan form 



Co 




(59) 
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We take the general form of Ci to be Ci = , I and impose the conditions 

\c d J 

and (1301). from which we obtain the constraints 



a + d = 2e (60) 

ad — be = 1. (61) 

The first constraint is solved by the parametrization a = e + u and d = e — u, which 
when inserted into the second equation, gives u^ + be = 0, the solution of which 
is b = uu and c = — 7. However the parameter u can be set to unity by a gauge 

transformation C,- — > UCjU^^ with U = \ ^ . Thus the final form of the 

matrices become 

One can verify the existence of both symmetries, with matrices 

X^( — ) ana ..(» i,). (03) 

such that 

x-^dx = Cj, n-^dn = cf. (64) 

The eigenvalues of the transfer matrix Ec = Cq ® Cq + Ci ® Ci^ are found to be 



ug , 1 



A5r=2, A^=2-n(7, A^ = 2 + ^ ± -^lOn^ + ^V- (65) 

Figure ^ shows the largest eigenvalue of the transfer matrix in the u — g plane. It 
is seen that the lines n = and g = Q are the crossover lines where the largest eigen- 
value changes and hence an MPS-QPT (MPS quantum phase transition) is expected 
to occur on these lines. 

Finally we can find the parent Hamiltonian of this model. The system of equations 
(J17p has a nontrivial solution for k = 2> and the solution space is spanned by the 
vectors: 



|ef) = |001) + |110) - |011) - |100) 

|ep) = (1 + ug){\mQ) + |111)) + (|001) + |011) + |100) + |110)) - 3(|010) + |101)) 

|ef) = 2(|000) -|111)) -3(|001) -|011) + |100) -|110)) 
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Figure 3: (Color Online) The largest eigenvalue of transfer matrix for model III. There 
are two crossover lines for the largest magnitude eigenvalue. 



K 



C\ 



2(|010) - |101)) - (1 + ug){\001) - |011) + |100) - |110)) 



(66) 



The interesting point about these vectors is that all of them are invariant (mod- 
ulo a sign) under the parity and spin-flip transformations and hence the symmetric 
Hamiltonian can be written in the following form, with four free coupling constants: 



he 



4 

E 



C\/„Ci 



JH){e^ 



(67) 



The explicit expression of the total Hamiltonian can be obtained along the same 
lines as for model A and model B. We will not do it here. 

7 Discussion 

In the formalism of matrix product states, there is a large room for constructing 
states and parent Hamiltonians. What really constrains this freedom and guides us 
along a way which may lead to interesting states and Hamiltonians is consideration 
of symmetries. The other constraining elements is the dimension of matrices which 
we choose. In this article we have classified all such states which are constructed from 
two dimensional matrices and have two important symmetries, namely the spin flip 
symmetry and the parity symmetry. We have shown that there are three different 
models, two of which lead to parent Hamiltonians with nearest and next-nearest 
interaction (models A and C) while one of them lead to a Hamiltonian with nearest 
neighbor interaction (model B). Furthermore, by calculating the eigenvalues of the 
transfer matrix in each case and determining the points of crossover between the 
largest and the next to largest eigenvalues of this matrix, we have identified the points 
of possible MPS quantum phase transitions. While in many of the works which have 
been reported on model building in matrix product states, [8l [9l [T0 | [T9 | [20 | [2T | [22 | [23] . 
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rotation symmetry has been taken into account, a condition which highly restricts 
the form of matrices, in this article we have relaxed this continuous symmetry in 
order to find all models compatible with discrete symmetries in order exhaust all the 
possibilities with two dimensional matrices. Any model with these symmetries must 
be equivalent to one of the above three models. For example in [27], the following 
model was suggested 

which has spin-flip symmetry with X = I n ]• "^^^ Hamiltonian for this model 

is [27] 

N 
H = J2 2(9' - l)Tt<jUi - (1 + g)Vf + iq- l)Vf af+iaf+2- (69) 

It is easy to see that this model is equivalent to model A above. In fact the transfor- 
mation 

SA{g = l)S-' = 2A[ 

with S = \ . n 2 in which q = ^i^^, proves the equivalence. 

V s™ 2 -cos 2 y 

8 Acknowledgements 

I would like to thank V. Karimipour for valuable discussion and comments. 

References 

[1] M. A. Nielsen, and I. L. Chnaivg] Quantum computation and quantum informa- 
tfon,Cambridge University Press, Cambridge, 2000. 

[2] M. Fannes, B. Nachtergaele and R. F. Werner, Commun. Math. Phys. 144, 443 
(1992). 

[3] I. Affleck, T. Kennedy, E. H. Lieb, H. Tasaki, Commun.Math. Phys. 115, 477 
(1988); 

[4] A. Klumper, A. Schadschneider and J. Zittartz, J. Phys. A (1991) L293; Z. Phys. 
B, 87 (1992) 281; Europhys. Lett., 24 (1993) 293. 

[5] H. Niggemann, and J. Zittartz, J. Phys. A: Math. Gen. 31, 9819-9828 (1998). 

[6] E. Bartel, A. Schadschneider and J. Zittartz, Eur. Phys. Jour. B, 31, 2, 209-216 

(2003). 



14 



[7] M. A. Ahrens, A. Schadschneider, and J. Zittartz, Europhys. Lett. 59 6, 889 
(2002). 

[8] A. K. Kolezhuk and H. J. Mikeska, Phys. Rev. Lett. 80, 2709 (1998); Int. J. Mod. 
Phys. B, 12, 2325-2348 (1998). 

[9] A. K. Kolezhuk, H. J. Mikeska, and Shoji Yamamoto, Phys. Rev. B 55, R3336 - 
R3339 (1997). 

[10] H. Niggemann, A. Kliimper, and J. Zittartz, Eur. Phys. J. B 13, 15 (2000). 

[11] J. M. Roman, G. Sierra, J. Dukelsky, and M. A. Martn-Delgado, J. Phys. A: 
Mathematical and General, 31, 48, 9729-9759(1998). 

[12] H. Niggemann, A. Kliimper, and J. Zittartz, Z. Phys. B 104, 103 (1997). 

[13] A. Kliimper, S. Matveenko, and J. Zittartz, Z. Phys. B 96, 401 (1995). 

[14] F. Verstraete, D. Porras, J. L Cirac, Phys. Rev. Lett. 93, 227205 (2004). 

[15] F. Verstraete, J. J. Garcia-RipoU, and J. I. Cirac, Phys. Rev. Lett. 93, 207204 
(2004). 

[16] F. Verstraete, J. I. Cirac, Phys. Rev. B 73, 094423 (2006); T. J. Osborne, Phys. 
Rev. Lett. 97, 157202 (2006); M. B. Hastings, Phys. Rev. B 73, 085115 (2006). 

[17] G. Vidal, Phys. Rev. Lett. 91, 147902 (2003). 

[18] F. Verstraete, M. A. Martin-Delgado, and J. L Cirac, Phys. Rev. Lett. 92, 087201 
(2004); F. Verstraete, M. Popp, and J. I. Cirac, Phys. Rev. Lett. 92, 027901 (2004); 
F. Verstraete, J. I. Cirac, J. L Latorre, E. Rico, and M. M. Wolf, Phys. Rev. Lett. 
94 140601(2005). 

[19] M. Asoudeh, V. Karimipour, and A. Sadrolashrafi, Phys. Rev. A 76, 012320 
(2007). 

[20] M. Asoudeh, V. Karimipour and A. Sadrolashrafi, Phys. Rev. B, 75, 224427 
(2007). 

[21] M. Asoudeh, V. Karimipour, and A. Sadrolashrafi, Phys. Rev. A 76, 012320 
(2007). 

[22] S. Alipour, V. Karimipour and L. Memarzadeh, Phys. Rev. A 75, 052322 (2007). 

[23] A. Trebedi and I. Bose, Phys. Rev. A 75, 042304(2007). 

[24] S. Alipour, S. Baghbanzadeh, and V. Karimipour, C Europhysics Letters (EPL), 
84 (2008) 67006. 

[25] V. Karimipour and L. Memarzadeh, Phys. Rev. B 77, 094416 (2008). 



15 



[26] S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cam- 
bridge, 1999). 

[27] M. M. Wolf, G. Ortiz, F. Verstraete and I. Cirac, Phys. Rev. Lett. 97, 110403 
(2006). 

[28] S. Alipour, V. Karimipour, and L. Memarzadeh, Eur. Phys. J. B 62, 159-169 
(2008). 



16 



